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Abstract 


Crop residues returned to the soil are important for the preservation of soil quality, health, and biodiversity, 
and they increase agriculture sustainability by recycling nutrients. Sugarcane is a bioenergy crop that produces 
huge amounts of straw (also known as trash) every year. In addition to straw, the ethanol industry also gener- 
ates large volumes of vinasse, a liquid residue of ethanol production, which is recycled in sugarcane fields as 
fertilizer. However, both straw and vinasse have an impact on N2O fluxes from the soil. Nitrous oxide is a 
greenhouse gas that is a primary concern in biofuel sustainability. Because bacteria and archaea are the main 
drivers of N redox processes in soil, in this study we propose the identification of taxa related with NO fluxes 
by combining functional responses (N2O release) and the abundance of these microorganisms in soil. Using a 
large-scale in situ experiment with ten treatments, an intensive gas monitoring approach, high-throughput 
sequencing of soil microbial 165 rRNA gene and powerful statistical methods, we identified microbes related to 
N20 fluxes in soil with sugarcane crops. In addition to the classical denitrifiers, we identified taxa within the 
phylum Firmicutes and mostly uncharacterized taxa recently described as important drivers of NxO consump- 
tion. Treatments with straw and vinasse also allowed the identification of taxa with potential biotechnological 
properties that might improve the sustainability of bioethanol by increasing C yields and improving N efficiency 


in sugarcane fields. 
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Introduction 


Bioethanol from sugarcane is becoming an increasingly 
important alternative energy source worldwide. The 
sugarcane planted area in Brazil is around 9 million-ha 
(CONAB, 2013) and 26 million-ha worldwide (Kao 
et al., 2014). The environmental benefits of replacing fos- 
sil fuels by ethanol from sugarcane can be only 
achieved if management practices are applied which 
lead to a minimum of greenhouse gases (GHG) losses to 
the atmosphere. Approximately 45% of the total amount 
of methane (CH4) and nitrous oxide (N20) emitted by 
the cultivation of sugarcane are derived from straw 
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burning (Macedo et al., 2008); this practice also releases 
charcoal particles, causing health problems for humans 
and animals. Recently, legislation was passed that 
restricts burning, which is being phased out. The 
unburned cane (green cane) leaves a thick mulch of 
plant material after harvest (trash or straw) ranging 
from 10 to 20 t ha~' (dry matter) consisting of leaves 
and tops. The conversion from burning to green man- 
agement of sugarcane will impact the biogeochemical 
cycling of carbon and nitrogen in the plant-soil system. 
Straw preservation affects the entire production process 
of sugarcane, influencing yields, fertilizer management 
and application, soil erosion, soil organic matter dynam- 
ics, and GHG emissions. 

In addition to plant straw, a second by-product of 
the sugarcane ethanol industry is vinasse (Mutton et al., 
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2010). Each liter of ethanol generates 10-15 1 of vinasse, 
which is a liquid organic residue that has a high poten- 
tial for polluting water streams. For this reason, there is 
legislation regulating the destination of vinasse. Cur- 
rently, Brazil produces approximately 337 billion liters 
of vinasse that is almost entirely recycled as fertilizer in 
sugarcane fields (Mutton et al., 2010). Vinasse has a high 
dissolved organic carbon (DOC) content resulting in 
high biochemical oxygen demand (BOD) and contains 
several nutrients; it is particularly rich in potassium and 
sulfur and contains considerable amounts of phospho- 
rus, nitrogen, calcium, magnesium, and micronutrients. 
Although vinasse use is of relevance for nutrient cycling 
in sugarcane fields its application contributes to GHG 
emissions especially when applied with chemical fertil- 
izers and on the straw (Carmo et al., 2013). 

The major sources of atmospheric NO are the nitrifi- 
cation and denitrification processes (Bouwman et al., 
2013). Globally, approximately 10.9 Tg N-N2O is esti- 
mated to be lost from land surfaces every year (Gallo- 
way et al., 2004). Although the metabolic pathways that 
result in N2O production during nitrification are not 
entirely clear, some authors suggest that this process 
might represent an even larger source of N2O from agri- 
cultural fields than denitrification (Khalil et al., 2004; 
Toyoda et al., 2011). Moreover, nitrification supplies oxi- 
dized N compounds that can be used as electron accep- 
tors, subsequently resulting in MNO release. 
Denitrification can also be a detoxification pathway as a 
result of the accumulation of NO, in the system (Zhu 
et al., 2013). Nitrogen oxidation regulates nutrient avail- 
ability for plant and the element mobility on the system, 
making the nitrification process of agronomical and 
environmental interest. Nitrification is a microbial- 
restricted niche (Hayatsu ef al., 2008) while denitrifica- 
tion can be performed by a broad diversity of microor- 
ganisms, especially prokaryotes (Philippot et al., 2007). 
Denitrification has being described as a modular path- 
way in which different cells may drive different steps to 
the complete reduction of nitrate to Nz being N2O an 
intermediate product (Zumft, 1997). Complementarily, 
the diversity of potential N,O reducers goes beyond 
sensu stricto denitrifiers (Jones ef al., 2013), which makes 
it difficult to identify the main players in N2O produc- 
tion and consumption. Identifying these players will 
improve our understanding of these processes in the 
field and promote more sustainable fertilizer use. 

While experiments on micro and mesocosm scales 
have been performed to identify microbial drivers of 
the N cycle and specifically linked to NO fluxes from 
soils (e.g. Ishii et al., 2011; Jung et al., 2014), field condi- 
tions cannot be reproduced. Surely such micro and mes- 
ocosm studies give important insights into the microbes 
involved in N cycle and specific processes. However, 


the dynamic system of soil/plants in addition to 
weather instability results in a complex microbial com- 
munity that goes beyond cultivable microorganisms. 
Thus, in this study by combining intensive in situ NO 
monitoring, high-throughput sequencing of soil micro- 
bial 16S rRNA gene and powerful statistical methods 
we have identified potential microbial taxa that might 
be explored to improve the sustainability of sugarcane 
as a bioenergy crop as a result of their biotechnological 
potential and/or their role in nutrient cycling. 


Materials and methods 


Site description and experimental design 


The experimental site is located in Piracicaba, Sao Paulo state, 
Brazil (22°41'19.34"S; 47°38’41.97"W). More than 50% of the 
sugarcane of Brazil is grown in this state (Conab, 2013). The 
site’s altitude is 575 m, the mean annual temperature is 21°C, 
and the mean annual precipitation is 1390 mm. The climate is 
defined as humid tropical with a rainy summer and a dry win- 
ter and is considered to be a Cwa type according to the 
Koppen classification system (Critchfield, 1960). The soil is clas- 
sified as Haplic Ferralsol (FAO, 2006); fertility parameters are 
listed in Table S1. Prior to the experiment, the sugarcane crop 
of the plant cycle was harvested without burning; thus, the 
field was in the first ratoon stage. After harvesting, the site was 
subdivided into two subsites. The straw was left onto soil of 
subsite 1 (10 t ha~') and removed mechanically of subsite 2. 
The removal of straw simulated the scenario that occurs when 
the sugarcane straw is used as raw material for biofuel produc- 
tion. 

Ten different treatments based on combination of organic, 
mineral fertilizers and straw were applied to the experimental 
site. The plots (n = 4; 8 x 9m) were randomly distributed 
and separated from each other by 2 m borders. All treatments 
received single superphosphate as the phosphorus source 
(17 kg ha! of P). The source of potassium was potassium 
chloride (100 kg ha~' of K) or vinasse (normal or concen- 
trated). The nitrogen source was ammonium nitrate applied 
at a rate of 100 kg N ha™'. The treatments without straw are 
described as the following: Ctrl (without N) — mineral fertil- 
izer containing P and K; MN (mineral N) — mineral fertilizer 
containing N, P, and K; MN_V (mineral N plus vinasse) — 
vinasse as K source (1.10° 1 ha~') and mineral fertilizer con- 


taining N and P; V (vinasse) — vinasse as K source 
(1.10° 1 ha~!) and mineral fertilizer containing P; CV (concen- 
trated vinasse) -— concentrated vinasse as K source 


(1.10 1 ha~!) and mineral fertilizer containing P. The same 
treatments were used with straw (Ctrl S, MN _ S; MN V S; 
V_S; and CV_S). Fertilization was carried out on 19 Novem- 
ber, 2012. Vinasse was applied to the total area, while the 
mineral fertilizers and concentrated vinasse were applied in 
bands parallel to the crop line as is usually performed in 
commercial areas. There was no treatment with CV and N 
because they could be applied to different sides of the plant; 
thus, the combination of C and N that results in high N20 
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emissions was avoided. The composition of the vinasse and 
concentrated vinasse is presented in Table $2. 


Greenhouse gas flux measurement and emission 
estimation 


Soil GHG fluxes were measured using the chamber-based 
method (Varner et al., 2003). In each plot, a PVC chamber 
(30 cm in diameter) was placed on the soil in which it was par- 
tially inserted (3 cm deep). After closing the chambers, 60 ml 
samples were collected at time points 1, 10, 20, and 30 min 
using syringes and stored under pressure in 20-ml-evacuated 
penicillin flasks sealed with gas-impermeable butyl-rubber 
septa (Bellco Glass 2048). All samplings were performed 
between 8:00 and 12:00 a.m. The first 31 samplings were col- 
lected within 60 days after fertilization, followed by 4 weekly 
samplings and finally fortnightly intervals until the crop was 
harvested. Each gas chamber flux was calculated from slope 
regression between the gas concentration and collection time 
according to Carmo et al. (2013). Measurements of atmospheric 
pressure, chamber height, and air temperature were taken dur- 
ing gas sampling to determine the air chamber volume and to 
calculate GHGs emissions. 

The gas samples were analyzed in a GC-2014 model gas 
chromatograph with electron capture for N2O (Shimadzu, 
Kyoto, Japan). Soil temperature was measured at a depth of 
10 cm to assist in the interpretation of the results. All of the 
chambers were installed at the fertilized sugarcane line posi- 
tion. GHG emissions were calculated by linear interpolation of 
fluxes between sampling events (Allen et al., 2010). We esti- 
mated that the fertilized bands accounted for 20% of the total 
experimental area; the remaining 80% of the area consisted of 
the spaces between the rows. Fertilizer emission factors (EF) 
were calculated according to Allen et al. (2010) and Carmo et al. 
(2013). 


Soil sampling and soil chemical analysis 


Soil sampling and gas collection initiated 1 day after fertiliza- 
tion. After collecting the gases, a composite soil sample 0- 
10 cm deep from three points was collected from each plot at 
the fertilizer band position. The samples were stored at —20°C. 
Soil moisture was determined by gravimetry, and mineral N 
content was determined by colorimetry of the soil extracts (KCl 
2 mol 17') using flow injection analysis (FIAlab 2500) based on 
the methods proposed by Kamphake et al. (1967) for NO; and 
Krom (1980) for NH. 


DNA extraction, 16S rRNA gene amplification and 
sequencing 


DNA was extracted from 250 mg of soil using the PowerSoil 
DNA Isolation Kit (Mobio Laboratories, Carlsbad, CA, USA) 
according to the manufacturer’s protocol. We performed DNA 
extraction in samples from eight time points of each treatment 
corresponding to 240 samples (10 treatments x 8 time 
points x 3 replicates). Archaeal and bacterial community struc- 


tures were assessed by sequencing the V4 region of the 16S 
rRNA gene using primers 515F (GIGCCAGCMGCCGCGG- 
TAA) and 806R (TAATCTWTGGGVHCATCAGG) (Caporaso 
et al., 2010) and the LIB-L kit for unidirectional sequencing. 
Duplicate 50 ul PCRs were performed as follows: 5 ul Roche 
10 x PCR buffer with MgCh, 2 U of Roche Tag DNA polymer- 
ase, 2 ul of dNTP (10 mM), 1 pl each primer (5 pmol l~’), and 
1 pl of DNA template. Thermocycling conditions consisted of 
an initial denaturation at 94 °C for 5 min, 30 cycles at 95 °C for 
30 s, 53 °C for 1 min, 72 °C for 1 min, and an final extension at 
72 °C for 10 min. Technical replicates were pooled and purified 
using the Qiagen PCR purification kit. The 240 samples were 
equimolar mixed and then sequenced (Macrogen Inc. Com- 
pany, South Korea) on a Roche 454 automated sequencer using 
the GS-FLX Titanium system (454 Life Sciences, Brandford, CT, 
USA). 


Sequence data analysis 


Sequence data were processed using the workflow of the NG6 
system (Mariette et al., 2012), installed on a local server, which 
system depends on Mothur version 1.31.1 (Schloss et al., 2009). 
The flowgrams were demultiplexed and sizes between 350 and 
750 flows were selected. The flowgrams were corrected using 
the shhh.flows command, which is the Mothur implementation 
of the original PyroNoise algorithm (Quince et al., 2011). Next, 
the results of the different sff files were combined for further 
analysis. The merged sequences were aligned to the bacterial 
and archaeal reference alignment provided on the Mothur web- 
site (http: //www.mothur.org/wiki/Silva_reference_alignment) 
based on the SILVA 102 release of the SSURef database (Quast 
et al., 2013). Only reads that fell into region 13 862-22 580 of 
the reference alignment were kept. The pre.cluster command 
was used to reduce typical sequence errors due to sequencing 
PCR products from high-diversity DNA samples. This com- 
mand assigns sequences that are within two mismatches to the 
most abundant sequence (Huse etal., 2010). Chimeric 
sequences were identified and removed using the chi- 
mera.uchime command (Edgar et al., 2011). Operational taxo- 
nomic units (OTUs) were formed at a maximum distance of 
0.03 using the dist.seqs command and average neighbor clus- 
tering. All sequences were taxonomically classified using the 
Mothur implementation of the RDP classifier (Wang et al., 
2007) using the training set (version 9) provided on the Mothur 
website — (http:/ /www.mothur.org/wiki/RDP_reference_files) 
and a bootstrap cutoff of 80%. OTUs were defined at an iden- 
tity cutoff of 97%. We are aware that the Mothur developers 
state this training set as of poor quality and it is small in size. 
However, RDP taxonomy is one of the only sets that take sub- 
groups of Acidobacteria into account. For each OTU, consensus 
taxonomy was determined using the classify.otu. command. 
Representative sequences for each OTU were re-aligned to the 
Silva reference alignment and a neighbor joining tree was cre- 
ated using the clear-cut program (Sheneman et al., 2006). Taxo- 
nomic classification and OTU clustering data were combined 
into the BIOM format (McDonald et al., 2012) for further down- 
stream statistical analysis with the Phyloseq (McMurdie & 
Holmes, 2013) R package. 
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The raw 454 pyrosequencing data of the 16S rRNA are avail- 
able at the European Nucleotide Archive (ENA) (https:// 
www.ebi.ac.uk/ena/) under the study Accession Number 
PRJEB8973. 


Statistical analyses 


Differential analyses of the count data (Anders & Huber, 2010) 
were applied to verify the effect of different management prac- 
tices on the abundance of operational taxonomic units (OTU). 
Time points were not considered for these analyses. The con- 
trol (Ctrl) treatment was used as the reference, and the negative 
binomial generalized linear model was applied using the 
DESeq2 package (Love et al., 2014) in R. OTU fold changes 
were considered statically significant using P < 0.01 as the 
criterion. 

Generalized linear mixed models (Bolker et al., 2009) were 
used to test the effect of different management practices, soil 
parameters, and soil taxa abundance on N30 fluxes. The sam- 
pling event was used as the random factor. First, we fitted a 
global model using the ‘Ime4’ R package version 1.1-7 (Bates 
et al., 2014). Model selection was performed using the ‘MuMIn’ 
R package version 1.10.5 (Barton, 2014) with the Akaike infor- 
mation criterion (AIC) for ranking. To verify the effect of OTUs 
on NO fluxes, we performed screening steps to reduce the 
number of parameters before applying the ‘dredge’ function, 
which imposes limited number of parameters. From each glo- 
bal model, we consecutively fitted nested models selecting the 
OTUs with P values <0.05, <0.1, and <0.2, until we obtained 
approximately 15 OTUs. Thus, the ‘dredge’ function was 
applied to select the best fitted model based on all parameter 
combinations. The AIC was used as a criterion for model selec- 
tion. In general, P values in generalized linear mixed models 
should be only used as guides to test nested models (Bolker 
et al., 2009). Thus, once the best model for each management 
practice was selected, and all OTUs were considered relevant 
for explaining N20 fluxes from soil even at P > 0.05. 


Results 


GHG fluxes and emissions 


Taking into account management practices and soil 
physicochemical factors, we found by the best fitted 
model that organic C addition (P < 0.001), N fertiliza- 
tion (P < 0.001), straw (P = 0.016) and moisture 
(P = 0.004) increased the fluxes of N2O from the soil. In 
contrast, NO; reduced (P = 0.065) the fluxes of N20 
from the soil. Ammonium (Fig. S1) and soil temperature 
did not explain N2O fluxes according to the best fitted 
model. The intercept of the model highly correlated 
with soil moisture (17 = 0.949). The relationship of 
nitrate with N20 production is well known. In our 
study, the models showed an association but not a 
cause-effect relationship because straw and organic fer- 
tilization increased the N2O fluxes but also reduced soil 
NO, (e.g., by immobilization; Fig. S2). Thus, N fertilization 


by itself (ie., available N and not the total N amount) bet- 
ter explained the N20 fluxes than soil inorganic N. The 
parameter estimation for fixed effects in global and best fit- 
ted models used to explain the N2O fluxes from soil con- 
sidering physicochemical parameters and management 
practices are presented in Table S3. Despite the high N2O 
peaks verified in treatments without straw (Fig. 1), annual 
N20 emissions were greater in the plots receiving treat- 
ment with straw (Fig. 2). We assumed that the effect of 
high moisture in the treatments with straw resulted in 
longer slopes with high N20 fluxes compared to treat- 
ments without straw (illustrated in Fig. 1). 

During the first 60 days of sampling, the cumulative 
N20 emissions reached more than 70% of the total emis- 
sions observed during the sugarcane cycle for all treat- 
ments (Fig. S3). At the 89th day, the cumulative 
emissions reached at least 85%. After 60 days, the soil 
inorganic N content reached a similar content to that 
before fertilization (Figs S1 and S2), which was expected 
as sugarcane plants were actively taking up nutrients 
from the soil. Therefore, we selected soil samples from 
the first 60 days to assess the microbial community 
related to N redox in soil. 

In accordance with the modeling showing that 
organic fertilization, mineral fertilization, and straw 
addition increased NO fluxes (Fig. 1), the calculated 
annual emissions were also affected by these manage- 
ment practices (Fig. 2). The emission factors (EF) 
showed that the amount of N from fertilizer released as 
N2O was equivalent to the following: (i) 210 g ha! 
when only mineral nitrogen was applied to the soil 
without straw; (ii) 2157 g ha! when mineral nitrogen 
was combined with vinasse; (iii) 1141 g ha! only with 
vinasse; and (iv) 503 g ha™™ when concentrated vinasse 
was used as fertilizer. In the treatments with straw, the 
N-N2O released from the fertilizers increased to 
1060 g ha! in the treatment containing mineral nitro- 
gen; 3380 g ha’! in the treatment with mineral nitrogen 
combined with vinasse; 1677 g ha! in the treatment 
with vinasse; and 688 g ha~' when concentrated vinasse 
was applied. 


Archaeal and Bacterial microbial community 


Overall sequencing of 165 rRNA gene marker yielded 
>1.3 million reads with an average length of 283 bp. 
After quality trimming, a total of 836 750 sequences 
were obtained from 240 samples. These sequences cov- 
ered 21 657 OTUs; a total of 7037 of these OTUs were 
represented by more than five reads. Differential analy- 
ses showed the effects of the management practices on 
423 OTUs. Time points were not discriminate for this 
purpose. From these, 32 belonged to Acidobacteria 
(Fig. 3c), 46 to Actinobacteria (Fig. 3f), 29 to Bacterioide- 
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N-N20 efflux (mg m? d7’) 


Soil moisture (%) 


Days after fertilization 


Fig. 1 Box-plot representation of nitrous oxide fluxes from soil with sugarcane crop after soil fertilization. Ctrl: no nitrogen addition; 
MN: mineral nitrogen fertilization; V: vinasse; CV: concentrated vinasse; S: straw left on the soil. *Daily measurements for samplings 
1-4; 3-4 times/week measurements for samplings 5-31; once a week measurements for samplings 32-35; fortnightly measurements 
for samplings 36-47 (just before the crop harvesting). The first sampling was 1 day after fertilization. 
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Fig. 2. Annual nitrous oxide emissions from soils with sugar- 
cane and respective emission factors (%). Bars are the standard 
errors. Reference lines are the calculated emissions from the 
control (Ctrl) treatments without straw (dotted line) and with 
straw (dashed line). Ctrl: no nitrogen addition; MN: mineral 
nitrogen fertilization; MN_V: mineral nitrogen plus vinasse fer- 
tilization; V: vinasse; CV: concentrated vinasse. Transversal 
lines: treatments with straw left on soil. 


tes (Fig. 3b), 24 to Chloroflexi (Fig. 3a), 58 to Firmicutes 
(Fig. 3d), 13 to Gemmatimonadetes (Fig. 3b), 3 to 
Nitrospirae (Fig. 3b), 8 to Planctomycetes (Fig. 3b), 124 


to Proteobacteria (Fig. 3e), 2 to Thaumarcheota (Fig. 3b), 
3 to Thermotogae (Fig. 3b), 5 to Verrucomicrobia 
(Fig. 3b), and 76 are unclassified bacteria or belong to 
other phyla. 

Interestingly, taxa classified as Chloroflexales were 
increased in abundance only in the treatments in which 
the straw was left on the soil. All OTUs in this order 
belonged to the genus Roseiflexus. In contrast, the differ- 
ential analysis showed that the members of Ktedobacte- 
rales were underrepresented in soils with straw. There 
was a significant effect of straw in reducing the abun- 
dance of 12 OTUs from the Ktedobacterales order, while 
in the plots without straw, only 4 OTUs were underrep- 
resented. These OTUs were underrepresented 23 times 
with straw and only five times without straw (Fig. 3a), 
suggesting an inhibitory effect of straw on Ktedobacte- 
rales members. 

The taxa recognized as ammonia oxidizers and 
affected by the treatments with straw and mineral N 
were from Thaumarcheota genus Nitrososphaera (OTU- 
123, OTU-333, Fig. 3b) and bacterial genus Nitrosospira 
(OTU-487, Fig. 3e, Fig. S4) (Nitrosomonadales order). 
Nitrososphaera OTUs were overrepresented only in treat- 
ments with straw, while Nitrosospira sp. OTU increased 
abundance in treatments with the addition of mineral 
N. The thousand most archaeal OTUs belong to the 
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Fig. 3 Effect of different treatments on OTUs belonging to Chloroflexi (a), Acidobacteria (c), Firmicutes (d), Proteobacteria (e), Actin- 
obateria (f), or other phyla (b). Numbers represent the OTU identification. Color groups at Family (a, d), Order (b, e, f), or Genus (c) 
level. Circle size indicates the fold change of the respective OTU when compared to the control treatment (Ctrl). Overrepresented 
OTUs are in the upper part of the plot and underrepresented OTU are in the below part of the plot. Ctrl_S: no N with straw; MN_S: 
mineral nitrogen with straw; MN_V_S: mineral nitrogen plus vinasse with straw; V_S: vinasse with straw; CV_S: concentrated vinasse 
with straw; CV: concentrated vinasse; V: vinasse; MN_V: mineral nitrogen plus vinasse; MN: mineral nitrogen. 


genus Nitrososphaera. Within the entire 16S rRNA data- 
set, there was a clear dominance of the OTU 487 classi- 
fied as ammonia oxidizer bacteria, which representative 
sequence presents 99% of identity with N. multiformis. 
The other member of the Nitrosomonadales order (OTU 
2796, Gallionellaceae) was low abundant. Among the 
nitrite oxidizers, only Nitrospira spp. (OTU-82, OTU-190, 
OTU-481, Fig. 3b) showed overrepresentation in treat- 
ments with straw and vinasse. The OTU-136 (Fig. 3d) 
classified as Ammoniphilus sp. also related to nitrogen 
cycle increased abundance in treatment with the addi- 
tion of vinasse. Ammoniphilus is obligatory oxalotrophic, 


haloalkalitolerant bacteria and requires high ammonium 
concentration for growth (Zaitsev et al., 1998). 

Vinasse organic fertilization resulted in an increase in 
37 Firmicutes taxa; however, the Lactobacillales order 
was the most responsive in terms of fold change 
(Fig. 3d). The Lactobacillales OTUs were also highly 
abundant based on the model intercept analysis (Fig. 
S5). Many OTUs (24) within Clostridiales were stimu- 
lated by the different treatments. However, 18 of these 
OTUs had low-positive fold changes with treatments 
including concentrated vinasse and straw left on the soil 
(Fig. 3d). Shifts in members of family Leuconostocaceae 
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only occurred with concentrated vinasse treatment, 
while shifts in members of family Alicyclobacillaceae 
were observed with normal vinasse treatment. OTU-215 
and OTU-818 belonging to genera Megasphaera and Dial- 
ister, respectively, were shifted in all treatments with 
normal vinasse (Fig. 3d). 

Regarding the phylum Actinobacteria, members of 
class Thermoleophilia were clearly overrepresented in 
treatments including the presence of straw. Order 
Micrococales from the class Actinobacteria was overrepre- 
sented in soils with the absence of straw, while the OTU 
belonging to the order Micronomosporales was overrepre- 
sented in soils with straw combined with mineral nitro- 
gen fertilization (Fig. 3f). 

In general, straw on the soil and the input of nutrients 
(organic or mineral) reduced the abundance of OTUs 
within phylum Acidobacteria. However, no pattern 
among specific Acidobacteria OTUs and fertilization 
practices was observed. The presence of straw resulted 
in 36 negative fold changes within 19 different OTUs; in 
treatment without straw, there were only 12 negative 
fold changes within 10 OTUs. Only 9 OTUs belonging 
to phylum Acidobacteria slightly increased their abun- 
dance in soil with straw (Fig. 3c). 

Overall, Proteobacteria was favored by organic fertil- 
ization and inhibited by the absence of straw (Fig. 3e). 
The treatments without straw resulted in 38 negative 
fold changes within 25 OTUs, while treatment with 
straw resulted in 19 negative fold changes within 14 
OTUs (Fig. 3e). In treatments with organic fertilization, 
Alphaproteobacteria OTU-52, OTU-71, and OTU-317 that 
were classified as Novosphingobium, Gluconacetobacter 
and Acetobacter, respectively, were highly overrepre- 
sented (Fig. 3e). In addition to organic fertilization with 
concentrated vinasse, members of order Rhodocyclales 
(Betaproteobacteria) needed soil covered with straw to be 
overrepresented. The ability of members of order Rhodo- 
cyclales to use complex organic sources and perform 
important steps of the nitrogen cycle, such as denitrifi- 
cation during waste treatment, has been previously 
described (Hesselsoe et al., 2009). The class Deltaproteo- 
bacteria is another group that was overrepresented in 
soil with straw. Additionally, the members of this group 
were underrepresented in soil without straw (Fig. 3e), 
indicating a condition of dependency on straw. Interest- 
ingly, there were different effects at the OTU level 
within the Anaeromyxobacter spp. (Deltaproteobacteria) in 
the different treatments. From this taxon, OTU-50 and 
OTU-434 were overrepresented in soil with straw, OTU- 
100 and OTU-248 were inhibited with the input of min- 
eral nitrogen, and OTU-289 had a positive fold change 
with the addition of concentrated vinasse. Among the 
members of Gammaproteobacteria, OTU-76 and OTU-120 
belonging to genera Acetinobacter and Dyella, respec- 


tively, were highly abundant in the treatments without 
straw when compared with the treatments with straw. 
Moreover, OTU-120, assigned to Dyella sp., was 
increased in soil with mineral nitrogen treatment with- 
out straw. 

OTU-121 and OTU-224 of the phylum Gemmatimona- 
detes were overrepresented in treatments with straw 
(Fig. 3b). Within the phylum Verrucomicrobia, OTUs 
belonging to the order Chthoniobacterales presented nega- 
tive fold changes in all treatments (Fig. 3b). One of these 
OTUs (OTU-398, classified as Candidatus Xiphinematobact- 
er) was underrepresented only in the treatments with 
straw. Members of this genus are obligate endosymbionts 
of nematodes of commercial importance (Vandekerckho- 
ve et al., 2000), suggesting that straw might reduce the 
amount of this parasite in soil. Therefore, this taxon might 
be considered as a bioindicator of soil quality. 


Microbial taxa to explain N20 fluxes 


After identifying the management practices (i.e., straw, 
mineral nitrogen fertilization, and organic fertilization) 
that increased N20 fluxes, we statistically checked what 
OTUs explain gas fluxes under the respective condi- 
tions. As a threshold, only OTUs that were represented 
at least twice in treatments with the same management 
practice were used to compose the respective model. 
These included were 64 OTUs related to straw, 24 OTUs 
related to mineral nitrogen, and 64 OTUs related to 
organic fertilization. In contrast to the differential analy- 
sis, we used frequency-normalized data to troubleshoot 
convergence failures. The global models that explain 
N20 fluxes for the different management practices using 
microbial taxa narrowed the field (P < 0.05) to 19 OTUs 
for straw, six for mineral nitrogen, and 14 for organic 
fertilization (Table 1). The best and global models fitted 
with microbial OTUs under different managements to 
explain N20 are presented in Table 1. The most repre- 
sented orders in these models were Burkholderiales (6 
OTUs), Bacillales (5 OTUs), Lactobacillales (4 OTUs), 
Myxococcales (3 OTUs), and Xanthomonadales (8 OTUs). 
All nested models generated to obtain approximately 15 
OTUs prior to application of the ‘dredge’ function (see 
Materials and methods for details) as well as the global 
and best models are presented in a supplemental file. 
The relative abundances of OTUs presented in Table 1 
are illustrated for Firmicutes (Fig. S6), Proteobacteria 
(Fig. S7) and other phyla (Fig. S8). There was an effect 
of vinasse on Lactobacillus, concentrated vinasse on Leu- 
conostoc and straw on Rummeliibacillus OTUs, signifi- 
cantly explaining the N,O fluxes from soil (Fig. S6). 
Furthermore, the effect of straw and vinasse on relative 
OTU abundance within Anaeromyxobacter spp. was 
prominent (Fig. S7). The effects of straw on the OTUs 
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Table 1 OTUs that explain N,O fluxes in soils under different sugarcane management practices (S = straw; MN = mineral nitrogen 
fertilization; V = vinasse fertilization) by the global models and by the best models 


Global model§ Best model Taxonomy 
OTU S MN V S MN V Phylum 
103 a Acidobacteria 
5 EAE % Actinobacteria 
231 * Actinobacteria 
429 * + Bacteroidetes 
503 ee Chloroflexi 
23 ii WE * BRE Firmicures 
61 ae F wE nae Firmicures 
436 bd Firmicures 
179 S = Firmicutes 
56 = Firmicutes 
136 AE ee * A A A Firmicutes 
144 * Firmicutes 
189 id *** Firmicutes 
366 ii PPE Firmicutes 
124 = i Firmicutes 
215 Hr Firmicutes 
121 me +*+ — Gemmatimonadetes 
224 ij Gemmatimonadetes 
34 e 7 Proteobacteria 
40 = Proteobacteria 
74 = + Proteobacteria 
88 * + f Proteobacteria 
341 * Proteobacteria 
826 *** Proteobacteria 
101 a Proteobacteria 
50 i = Proteobacteria 
289 a k Proteobacteria 
434 + Proteobacteria 
317 * Proteobacteria 
641 kil ** Proteobacteria 
131 ae Proteobacteria 
120 * bid Proteobacteria 
373 sas sa Proteobacteria 
469 t Proteobacteria 
333 * Thaumarchaeota 
172 t Unclassified 
206 * ae Unclassified 
242 ae . $ + t Unclassified 


Order Family Genus 
Acidobacteria Gp6 
Actinomycetales Micrococcaceae Arthrobacter 
Solirubrobacterales Unclassified 
Sphingobacteriales Sphingobacteriaceae Sphingobacterium 
Chloroflexales Chloroflexaceae Roseiflexus 
Lactobacillales Lactobacillaceae Lactobacillus 
Lactobacillales Lactobacillaceae Lactobacillus 
Lactobacillales Lactobacillaceae Lactobacillus 
Lactobacillales Leuconostocaceae Leuconostoc 
Bacillales Bacilaceae Bacillus 
Bacillales Paenibacillaceae Ammoniphilus 
Bacillales Paenibacillaceae Rummeliibacillus 
Bacillales Bacilaceae Unclassified 
Bacillales Sporolactobacillaceae Sporolactobacillus 
Clostridiales Clostridiaceae Clostridium 
Selenomonadales Veillonellaceae Megasphaera 
Gemmatimonadales Gemmatimonadaceae Gemmatimonas 
Gemmatimonadales Gemmatimonadaceae Gemmatimonas 
Burkholderiales Comamonadaceae Unclassified 
Burkholderiales Burkholderiaceae Burkholderia 
Burkholderiales Oxalobacteraceae Unclassified 
Burkholderiales Burkholderiaceae Burkholderia 
Burkholderiales Burkholderiaceae Ralstonia 
Burkholderiales Unclassified 
Enterobacteriales Enterobacteriaceae Unclassified 
Myxococcales Cystobacteraceae Anaeromyxobacter 
Myxococcales Cystobacteraceae Anaeromyxobacter 
Myxococcales Cystobacteraceae Anaeromyxobacter 
Rhodospirillales Acetobacteraceae Acetobacter 
Sphingomonadales Sphingomonadaceae Sphingomonas 
Unclassified 
Xanthomonadales Xanthomonadaceae Dyella 
Xanthomonadales Xanthomonadaceae Arenimonas 
Xanthomonadales Sinobacteraceae JTB255 
Nitrososphaerales Nitrososphaeraceae Nitrososphaera 


§Only OTUs that presented P < 0.05 are shown. Complete list of OTUs are presented in a supplemental file. 
***p < 0.001; **P < 0.01; *P < 0.05; +P < 0.1; {not significant but is part of the best model selected by AIC. 


that explain N20 fluxes within the archaeal genus Nitro- 
sosphaera and bacterial genera Roseiflexus and Gemmati- 
monas are shown in Fig. S8. 


Discussion 


Based on our field experiment, we expect to guide fur- 
ther efforts for mining specific members in the microbial 
community for sustainable biofuel production. A prom- 


ising example for mining is our finding of the domi- 
nance of archaeal WNitrososphaera and bacterial 
Nitrosospira as ammonia oxidizers instead of Nitrosomon- 
as commonly found in soils and used as model for 
ammonia oxidizing bacteria. One of the practices to 
reduce N20 emissions in sugarcane is the use of nitrifi- 
cation inhibitors. However, the available inhibitors are 
expensive and some are applied in high amounts (up to 
5% of the N is applied as nitrification inhibitor com- 
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pounds) into soil (Soares et al., 2015). Furthermore, the 
available nitrification inhibitors were initially addressed 
to Nitrosomonas (Slangen & Kerkhoff, 1984) and some of 
them have different effect on Nitrososphaera and Nitroso- 
spira (Shen et al., 2013). Therefore, based on our results, 
the use and development of specific products for these 
taxa might be addressed. 

The genus Roseiflexus were present in all treatments, 
however, overrepresented in treatments with straw 
(Fig. 3a). Of the six Chloroflexi classes, Chloroflexia are 
able to perform photosynthesis (Gupta et al., 2013). 
Despite the fact that mixotrophy metabolism have been 
reported for Roseiflexus spp. (Klatt et al., 2013), it is diffi- 
cult to suggest that photosynthesis occurred below the 
straw layer due to the absence/low incidence of light. 
Until recently, no other morphological, physiological, 
biochemical, or molecular traits were known to be 
shared by different groups from phylum Chloroflexi 
(Gupta et al., 2013). Nevertheless, Dellas et al. (2013) 
have described, using Roseiflexus sp. as reference organ- 
ism, that the members of class Chloroflexia present a sin- 
gular mevalonate pathway. In other words, they present 
a singular pathway to synthesize the precursor of all 
isoprenoids, which are compounds of biotechnological 
interest for a broad range of industrial areas (Maury 
et al., 2005). Additionally, some results have shown that 
members of phylum Chloroflexi are promising cellulose 
degraders (Stott et al., 2009) and the utilization of plant 
polymers may be widespread across the Chloroflexi 
phylum (Hug et al., 2013). Many efforts have being 
addressed to obtain cellulases by screening and 
sequencing new organisms as they have potential to 
become the largest group of industrially used enzymes 
worldwide (Wilson, 2009). These results highlight Ros- 
eiflexus as an important target for optimization of the C 
yield per harvested area from crops used for biofuel 
production. 

Large amounts of N are exported from the sugarcane 
crop as compared to the N applied as fertilizer in field, 
what makes this crop one of the most efficient crops for 
bioenergy. One of the explanations is that sugarcane 
benefits from N>-fixing processes (Boddey et al., 2003; 
Urquiaga et al., 2012). Many bacterial isolates and bacte- 
ria consortia have been used to test this hypothesis 
without consistent results (Boddey et al., 2003; Martinez 
De Oliveira et al., 2006; Schultz et al., 2014). The reason 
might be that a great number of N>-fixing organisms 
exist in the field, both in soils and in plants, non-endo- 
phytic microorganisms contribute to N fixation (Fischer 
et al., 2012). Although nitrogenase activity has not yet 
been confirmed in Roseiflexus spp. strains, transcriptom- 
ic profiles similar to Oscillochloris trichoides have shown 
that Roseiflexus spp. fix atmospheric N and this activity 
is more intense during the night (Klatt et al., 2013). 


Thus, deeper studies concerning soil living Roseiflexus 
spp. and their functional traits would contribute to bet- 
ter understanding of N cycle in soils. 

Here, we showed that organic fertilization was the 
management practice that most affected the abundance 
of Firmicutes in soils (Fig. 3d). As proposed by Shar- 
min et al. (2013), Firmicutes-type bacteria play a domi- 
nant role in the sugarcane industry due to their 
physiological ability to ferment sugars. Firmicutes, in 
particular lactic acids producers, are historically recog- 
nized as contaminants by the ethanol industry because 
their metabolites effectively inhibit yeast ethanol pro- 
duction (Muthaiyan et al., 2011). Thus, because most of 
the OTUs of Lactobacillus and Leuconostoc were not 
observed in treatments without vinasse application and 
Lactobacillus is the dominant taxa in vinasse (Assis 
Costa et al., 2015), vinasse residues might be the main 
source of these taxa to the soil (Fig. S6). As verified in 
this work and by Carmo et al. (2013), vinasse increases 
N2O emissions from sugarcane fields. The main reason 
might be because vinasse residues are rich in labile C 
and have high biological oxygen demand (BOD) values 
(Carmo et al., 2013; Christofoletti et al., 2013); thus, car- 
bon and moisture are the key factors related with NO 
release. General fermenters and Lactobacillus are pres- 
ent in vinasse and when vinasse is applied to soil 
those microorganisms might be additional important 
factor for the N2O balance. In accordance with Philip- 
pot et al. (2011), we also propose that the increase in 
the proportion of these organisms in soils here studied 
might result in increases in NO release as truncated 
denitrification pathways are found in Lactobacillus spp. 
(Shapleigh, 2013). Four of the taxa that explained NO 
fluxes from soil in our study were probably added to 
soil by vinasse organic fertilization. These taxa are 
from the genera Lactobacillus and Leuconostoc and were 
not detected in treatments without vinasse application 
(Fig. S6). 

Another important observation is that nitric oxide 
(NO) is a natural product from amino acids reduction 
during fermentation processes as well from amino 
acids oxidation (Medinets et al., 2015). NO can also be 
used as an electron acceptor by microorganisms such 
as the nondenitrifier Anaeromyxobacter spp., which have 
the ability to reduce NO and NO (Sanford et al., 
2012). Thus, the question arises of whether the propor- 
tion of nitrous oxide reducers can improve the NO 
balance. Nishizawa et al. (2014) showed that the inocu- 
lation into the soil of No-generating denitrifiers (Azoar- 
cus sp., Niastella sp. and Burkholderia sp.) reduced the 
release of N20. Nevertheless, it is not adequate in 
fields with mineral fertilization because denitrifiers 
also tend to reduce the efficiency of the applied N. 
One promising approach is the overexpression of the 
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nosZ gene by denitrifiers to improve the N20 balance, 
for instance, when N is not dependent on artificially 
fixed sources (Itakura et al., 2013). In theory, microor- 
ganisms lacking genes related to NxO generation (NO 
and N,O) or those with a greater affinity for NxO than 
other N species could also be applied to fertilized 
fields (Miyahara et al., 2010). Sanford et al. (2012) pro- 
posed that the proportion of nitrous oxide reducers 
would improve the N20 balance between the soil and 
the atmosphere, with special attention paid to the non- 
denitrifier Anaeromixobacter spp. This genus was among 
the most abundant taxa identified in this study (OTUs- 
50, -100, -133, -150, -236, -289, -434, -764, -1051, -1435, 
and -1593; Figs S5 and S7). Microorganisms that use 
trace gas as nutrient source or electron acceptors 
require efficient enzymatic systems because the concen- 
tration of these gases is very low (Conrad, 1996). By 
consensus, the more specialized the function is, the 
more diversity is required as the rare (micro) organ- 
isms exert these functions as has been shown for meth- 
ane consumption (Pester et al., 2010). Thus, it is not 
clear whether trace gas consumption is a specialized 
niche or whether diversity is not required to maintain 
active specialized niches. Because Anaeromyxobacter 
spp. can use sources other than NxO as an electron 
acceptor (i.e., iron oxides; Hori et al., 2010), we assume 
that the reduction of trace gases is not a specialty of 
Anaeromyxobacter. According to the thermodynamic 
theory, oxygen>nitrogen>manganese>iron>sulfur oxi- 
des are reduced sequentially (Ponnamperuma, 1972; 
Thauer et al., 1977). Therefore, in theory when NxO 
are available in the environment, Anaeromyxobacter spp. 
tend to use these gases rather than other oxides. Thus, 
NxO reduction does not represent the niche of these 
microorganisms and the reduction of NxO does not 
require specialized microorganisms. OTUs classified as 
Anaeromyxobacter spp. were identified both by the glo- 
bal models and by the best fitted models explaining 
N30 fluxes from soil with straw and soil with mineral 
N fertilization (Table 1). To date, there are four strains 
of Anaeromyxobacter spp. with available genome 
sequences that are able to reduce NO and N30; how- 
ever, they do not possess the capacity to produce NO 
from NO, (Sanford et al., 2012). Our findings linking 
this taxon to N20 fluxes reinforce the hypothesis that 
Anaeromyxobacter spp. are important contributors to 
N:O consumption. However, in this present study, 
some OTUs of Anaeromyxobacter spp. were overrepre- 
sented, while others were underrepresented under the 
management practices that increase N?O emissions. 
These results indicate that the selection of groups 
related to N2O balance at different taxonomical levels 
should be considered. For instance, our results suggest 
that different OTUs might share the same ability to 


perform a specific function, but they might respond dif- 
ferently under conditions of stress (i.e., fertilizers). 

Most OTUs identified in this study seem to be related 
to processes involved in N2O release. For example, one 
OTU of Nitrososphaera significantly explained the fluxes 
of N2O on the global model in straw management 
(Table 1, Fig. 3b, Fig. S8). This archaeal genus is well 
recognized as an ammonia oxidizer (Spang et al. 2010) 
with the capacity for N2O production from ammonia 
oxidation (Jung et al., 2014; Stieglmeier et al., 2014). 

Jones et al. (2014) showed that 33% of the nosZ gene 
(encoding nitrous oxide reductase) sequences in diverse 
environmental samples were most closely related to 
Gemmatinoma ssp., while 18% of all clones grouped with 
nosZ sequences were from phyla Bacteroidetes and 
Chloroflexi. These groups are members of the poorly 
characterized Clade II of the nitrous oxide-reducing 
microbial community (Jones et al., 2014), and they were 
represented in our fitted models explaining N20 fluxes 
from sugarcane soil. However, due to their poor level of 
characterization and lack of experimental data tracing 
N2O production or consumption such as for Anaero- 
myxobacter spp., it is not possible to support a cause- 
effect relationship between these taxa and N20 fluxes. It 
is also not possible to infer if they might contribute to 
N20 production or consumption. In contrast, the well- 
characterized denitrifying Proteobacteria were identified 
within the OTUs explaining N2O fluxes. For instance, 
strains of Burkholderia and Dyella have been tested to 
reduce the amount of N20 released into the atmosphere 
after organic fertilization (Nishizawa et al., 2014). How- 
ever, they might increase the N losses from the soil to 
atmosphere through denitrification and consequently 
reduce synthetic fertilizer efficiency. The OTU classified 
as Ammoniphilus sp. significantly explained N2O fluxes 
in all fitted models (Table 1). Despite the fact that this 
genus is known to require high concentrations of 
ammonium and use oxalates exclusively as a C source, 
to date no metabolic pathways that could result in N20 
have been described for members of this genus. The 
available Ammoniphilus sp. type strains are not able to 
reduce NO, and are strictly aerobic (Zaitsev et al., 
1998). 

In an ideal context, a nitrogen fixer such as Azoarcus 
sp. that lacks the steps of denitrification prior to the pro- 
duction of N2O (Krause et al., 2006) might be used to 
improve the N balance in bioenergy crops. Azoarcus 
strains have been tested for potential reduction of NxO 
emissions from rice fields (Ishii et al., 2011; Nishizawa 
et al., 2014). For sugarcane, an analogous situation could 
arise with members of the family Acetobacteraceae. All of 
the known N-fixing organisms that are members of this 
family belong to the genera Acetobacter and Glu- 
conacetobacter (Dutta & Gachhui, 2006). The association 
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between these two genera and sugarcane was reported 
by James et al. (1994). These organisms are of special 
interest for sugarcane crop as they have the ability to fix 
N in aerobic conditions like well-drained soils. Addi- 
tionally, Acetobacter and Gluconacetobacter can use a 
broad range of organic sources, including alcohols pres- 
ent in the vinasse (James et al., 1994; Parnaudeau et al., 
2008). In our study, OTUs (OTU-71 and OTU-317, 
Fig. 3e) belonging to these two genera were strongly 
favored by vinasse organic fertilization and as showed 
by Assis Costa et al. (2015) they are also present in the 
ethanol production process. The same Acetobacter OTU 
favored by vinasse organic fertilization significantly 
explained N,O fluxes from soil with straw by the global 
fitted model. Thus, Acetobacter and Gluconacetobacter 
should be highlighted as the closest representatives to 
an ideal situation for the improvement of N-cycling in 
sugarcane, both by acting as N-fixing organisms and 
potential drivers of processes related to NO release. 
Furthermore, Acetobacter and Gluconacetobacter could be 
added to vinasse used as organic fertilizer in sugarcane 
fields, and thus, contributing to redefine vinasse from a 
source to a sink of Nə and NxO. However, it must be 
acknowledged that the survival of microrganisms inocu- 
lated in the soil system is a challenge. 

To our knowledge, this is the first study relating spe- 
cific microbial taxa to the respective functional response 
(N20 release) along the time and under different man- 
agement practices in tropical soil. Here, we show that 
the limitation of linking microbial taxa with functional 
traits is particularly challenging because the microbial 
responses are at different taxonomical levels such as at 
OTU and genus levels, as observed for Anaeromyxobacter 
and Roseiflexus, respectively. 

The effects of fertilization and straw on NO emis- 
sions have been previously reported (Carmo etal., 
2013), however, not taken into account the microbial 
community. Besides the common Proteobacteria model 
organisms for denitrification, we have found taxa 
recently described as potential drivers of N2O produc- 
tion and consumption (Jones et al., 2013). Additionally, 
we have identified taxa with potential biotechnological 
properties that might improve the sustainability of bio- 
ethanol by increasing C yields and improving N effi- 
ciency in sugarcane fields. 


Acknowledgements 


We thank AC Vitti, AM Yasuda, KS Lourenço, L Martinelli for 
technical assistance and W van der Werf, L Hemerik for data 
analyses instructions. This study was supported by grants from 
CAPES/NUFFIC (037/2012; BEX 0258/13-0); Sao Paulo 
Research Foundation-FAPESP (2012/50694-6) and Joint 
Research Projects in Biobased Economy (FAPESP/NWO 13/ 
50365-5); The National Council of Technological and Scientific 


Development (CNPq) (471886/2012-2). Publication 5868 of the 
Netherlands Institute of Ecology (NIOO-KNAW). 


References 


Allen DE, Kingston G, Rennenberg H, Dalal RC, Schmidt S (2010) Effect of nitrogen 
fertilizer management and waterlogging on nitrous oxide emission from subtropi- 
cal sugarcane soils. Agriculture Ecosystems & Environment, 136, 209-217. 

Anders S, Huber W (2010) Differential expression analysis for sequence count data. 
Genome Biology, 11, R106. 

Assis Costa OY, Souto BM, Tupinamba DD et al. (2015) Microbial diversity in sugar- 
cane ethanol production in a Brazilian distillery using a culture-independent 
method. Journal of Industrial Microbiology & Biotechnology, 42, 73-84. 

Barton K (2014) MuMIn: Multi-model inference. [Computer software manual]. Avail- 
able at: http://cran.rproject.org/web/packages/MuMIn/MuMIn.pdf (accessed 
04 February 2014). 

Bates D, Maechler M, Bolker B, Walker S (2014) Ime4: Linear mixed-effects models 
using Eigen and S4. [Computer software manual]. Available at: http://cran.r- 
project.org/web/ packages /Ime4/Ime4.pdf (accessed 03 February 2014). 

Boddey RM, Urquiaga S, Alves BJR, Reis V (2003) Endophytic nitrogen fixation in 
sugarcane: present knowledge and future applications. Plant and Soil, 252, 139- 
149. 

Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MHH, White J-SS 
(2009) Generalized linear mixed models: a practical guide for ecology and evolu- 
tion. Trends in Ecology & Evolution, 24, 127-135. 

Bouwman AF, Beusen AHW, Griffioen J et al. (2013) Global trends and uncertainties 
in terrestrial denitrification and NO emissions. Philisophical Transactions of the 
Royal Society of London, B Biological Sciences, 386, 20130112. 

Caporaso JG, Lauber CL, Walters WA ef al. (2010) Global patterns of 165 rRNA 
diversity at a depth of millions of sequences per sample. Proceedings of the 
National Academy of Sciences of the United States of America, 108, 4516-4522. 

Carmo JB, Filoso S, Zotelli LC et al. (2013) Infield greenhouse gas emissions from 
sugarcane soils in Brazil: effects from synthetic and organic fertilizer 
application and crop trash accumulation. Global Change Biology Bioenergy, 5, 
267-280. 

Christofoletti CA, Escher JP, Correia JE, Urbano Marinho JF, Fontanetti CS (2013) 
Sugarcane vinasse: environmental implications of its use. Waste Management, 33, 
2752-2761. 

Conab (2013) Acompanhamento da Safra Brasileira: Cana-de-açúcar, safra 2013/2014. 
Companhia Nacional de Abastecimento, Brasília, 19 p. 

Conrad R (1996) Soil microorganisms as controllers of atmospheric trace gases (Ho, 
CO, CH;,, OCS, N20, and NO). Microbiological Reviews, 60, 609-640. 

Critchfield HJ (1960) General climatology. Prentice-Hall, Englewood Cliffs. 465 p. 

Dellas N, Thomas ST, Manning G, Noel JP (2013) Discovery of a metabolic alterna- 
tive to the classical mevalonate pathway. eLife, 2, e0067. 

Dutta D, Gachhui R (2006) Novel nitrogen-fixing Acetobacter nitrogenifigens sp nov., 
isolated from Kombucha tea. International Journal of Systematic and Evolutionary 
Microbiology, 56, 1899-1903. 

Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R (2011) UCHIME improves sen- 
sitivity and speed of chimera detection. Bioinformatics, 27, 2194-2200. 

FAO (2006) World Reference Base for Soil Resources, 2006: A Framework for International 
Classification, Correlation and Communication. Food and Agriculture Organization 
of the United Nations, Rome. 

Fischer D, Pfitzner B, Schmid M et al. (2012) Molecular characterisation of the diazo- 
trophic bacterial community in uninoculated and inoculated field-grown sugar- 
cane (Saccharum sp.). Plant and Soil, 356, 83-99. 

Galloway JN, Dentener FJ, Capone DG et al. (2004) Nitrogen cycles: past, present 
and future. Biogeochemistry, 70, 153-226. 

Gupta RS, Chander P, George S (2013) Phylogenetic framework and molecular sig- 
natures for the class Chloroflexi and its different clades; proposal for division of 
the class Chloroflexi class. nov into the suborder Chloroflexineae subord. nov., con- 
sisting of the emended family Oscillochloridaceae and the family Chloroflexaceae 
fam. nov., and the suborder Roseiflexineae subord. nov., containing the family 
Roseiflexaceae fam. nov. Antonie Van Leeuwenhoek International Journal of General and 
Molecular Microbiology, 103, 99-119. 

Hayatsu M, Tago K, Saito M (2008) Various players in the nitrogen cycle: diversity 
and functions of the microorganisms involved in nitrification and denitrification. 
Soil Science and Plant Nutrition, 54, 33—45. 

Hesselsoe M, Fuereder S, Schloter M et al. (2009) Isotope array analysis of Rhodo- 
cyclales uncovers functional redundancy and versatility in an activated sludge. 
ISME Journal, 3, 1349-1364. 


© 2015 The Authors Global Change Biology Bioenergy Published by John Wiley & Sons Ltd., doi: 10.1111/gcbb.12284 


12 L.M. PITOMBO et al. 


Hori T, Mueller A, Igarashi Y, Conrad R, Friedrich MW (2010) Identification of iron- 
reducing microorganisms in anoxic rice paddy soil by C-13-acetate probing. ISME 
Journal, 4, 267-278. 

Hug LA, Castelle CJ, Wrighton KC et al. (2013) Community genomic analyses con- 
strain the distribution of metabolic traits across the Chloroflexi phylum and indi- 
cate roles in sediment carbon cycling. Microbiome, 1, 22-22. 

Huse SM, Welch DM, Morrison HG, Sogin ML (2010) Ironing out the wrinkles in the 
rare biosphere through improved OTU clustering. Environmental Microbiology, 12, 
1889-1898. 

Ishii S, Ohno H, Tsuboi M, Otsuka S, Senoo K (2011) Identification and isolation of 
active N20 reducers in rice paddy soil. ISME Journal, 5, 1936-1945. 

Itakura M, Uchida Y, Akiyama H et al. (2013) Mitigation of nitrous oxide emissions 
from soils by Bradyrhizobium japonicum inoculation. Nature Climate Change, 3, 208- 
212. 

James EK, Reis VM, Olivares FL, Baldani JI, Dobereiner J (1994) Infection of sugar- 
cane by the nitrogen-fixing bacterium Acetobacter diazotrophicus. Journal of Experi- 
mental Botany, 45, 757-766. 

Jones CM, Graf DRH, Bru D, Philippot L, Hallin S (2013) The unaccounted yet abun- 
dant nitrous oxide-reducing microbial community: a potential nitrous oxide sink. 
ISME Journal, 7, 417-426. 

Jones CM, Spor A, Brennan FP et al. (2014) Recently identified microbial guild medi- 
ates soil N20 sink capacity. Nature Climate Change, 4, 801-805. 

Jung M-Y, Well R, Min D et al. (2014) Isotopic signatures of N2O produced by 
ammonia-oxidizing archaea from soils. ISME Journal, 8, 1115-1125. 

Kamphake LJ, Hannah SA, Cohen JM (1967) Automated analysis for nitrate by 
hydrazine reduction. Water Research, 1, 205-216. 

Kao MCJ, Gesmann M, Gheri F (2014) FAOSTAT: A complementary package to the 
FAOSTAT database and the Statistical Yearbook of the Food and Agricultural 
Organization of the United Nations [Computer software manual]. Available at: 
http: / /cran.r-project.org / web/packages/FAOSTAT/index.html (accessed 02 
June 2014). 

Khalil K, Mary B, Renault P (2004) Nitrous oxide production by nitrification and 
denitrification in soil aggregates as affected by O-2 concentration. Soil Biology & 
Biochemistry, 36, 687-699. 

Klatt CG, Liu ZF, Ludwig M, Kuhl M, Jensen SI, Bryant DA, Ward DM (2013) Tem- 
poral metatranscriptomic patterning in phototrophic Chloroflexi inhabiting a 
microbial mat in a geothermal spring. ISME Journal, 7, 1775-1789. 

Krause A, Ramakumar A, Bartels D et al. (2006) Complete genome of the mutualis- 
tic, N-2-fixing grass endophyte Azoarcus sp strain BH72. Nature Biotechnology, 24, 
1385-1391. 

Krom MD (1980) Spectrophotometric determination of ammonia — a study of a mod- 
ified Berthelot reaction using salicylate and dichloroisocyanurate. Analyst, 105, 
305-316. 

Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dis- 
persion for RNA-Seq data with DESeq2. Genome Biology, 15, 550. 

Macedo IC, Seabra JEA, Silva JEAR (2008) Green house gases emissions in the pro- 
duction and use of ethanol from sugarcane in Brazil: the 2005/2006 averages and 
a prediction for 2020. Biomass and Bioenergy, 32, 582-595. 

Mariette J, Escudie F, Allias N, Salin G, Noirot C, Thomas S, Klopp C (2012) NG6: 
integrated next generation sequencing storage and processing environment. BMC 
Genomics, 13, 462. 

Martinez De Oliveira AL, Canuto EDL, Urquiaga S, Reis VM, Baldani JI (2006) Yield 
of micropropagated sugarcane varieties in different soil types following inocula- 
tion with diazotrophic bacteria. Plant and Soil, 284, 23-32. 

Maury J, Asadollahi MA, Moller K, Clark A, Nielsen J (2005) Microbial isoprenoid 
production: an example of green chemistry through metabolic engineering. Bio- 
technology for the Future, 100, 19-51. 

McDonald D, Clemente J, Kuczynski J et al. (2012) The Biological Observation Matrix 
(BIOM) format or: how I learned to stop worrying and love the ome-ome. Giga- 


Science, 1, 7. 

McMurdie PJ, Holmes S (2013) phyloseq: an R package for reproducible interactive 
analysis and graphics of microbiome census data. PLoS ONE, 8, e61217. 

Medinets S, Skiba U, Rennenberg H, Butterbach-Bahl K (2015) A review of soil NO 
transformation: associated processes and possible physiological significance on 
organisms. Soil Biology and Biochemistry, 80, 92-117. 

Miyahara M, Kim S-W, Fushinobu S et al. (2010) Potential of aerobic denitrification 
by Pseudomonas stutzeri TR2 to reduce nitrous oxide emissions from wastewater 
treatment plants. Applied and Environmental Microbiology, 76, 4619-4625. 

Muthaiyan A, Limayem A, Ricke SC (2011) Antimicrobial strategies for limiting bac- 
terial contaminants in fuel bioethanol fermentations. Progress in Energy and Com- 
bustion Science, 37, 351-370. 


Mutton MA, Rossetto R, Mutton MJR (2010) Agricultural use of stillage. In: Sugar- 
cane Bioethanol - R&D for Productivity and Sustainability (ed. Cortez LAB), pp. 423- 
440. Editora Edgard Blücher, Sao Paulo. 

Nishizawa T, Quan A, Kai A et al. (2014) Inoculation with N-2-generating denitrifier 
strains mitigates N2O emission from agricultural soil fertilized with poultry man- 
ure. Biology and Fertility of Soils, 50, 1001-1007. 

Parnaudeau V, Condom N, Oliver R, Cazevieille P, Recous S (2008) Vinasse organic 
matter quality and mineralization potential, as influenced by raw material, fer- 
mentation and concentration processes. Bioresource Technology, 99, 1553-1562. 

Pester M, Bittner N, Deevong P, Wagner M, Loy A (2010) A ‘rare biosphere’ micro- 
organism contributes to sulfate reduction in a peatland. ISME Journal, 4, 1591- 
1602. 

Philippot L, Hallin S, Schloter M (2007) Ecology of denitrifying prokaryotes in agri- 
cultural soil. Advances in Agronomy, 96, 249-305. 

Philippot L, Andert J, Jones CM, Bru D, Hallin S (2011) Importance of denitrifiers 
lacking the genes encoding the nitrous oxide reductase for N20 emissions from 
soil. Global Change Biology, 17, 1497-1504. 

Ponnamperuma FN (1972) The Chemistry of Submerged Soils. Advances in Agron- 
omy, 24, 29-96. 

Quast C, Pruesse E, Yilmaz P et al. (2013) The SILVA ribosomal RNA gene database 
project: improved data processing and web-based tools. Nucleic Acids Research, 41, 
590-596. 

Quince C, Lanzen A, Davenport RJ, Turnbaugh PJ (2011) Removing noise from py- 


rosequenced amplicons. BMC Bioinformatics, 12, 38. 

Sanford RA, Wagner DD, Wu Q et al. (2012) Unexpected nondenitrifier nitrous oxide 
reductase gene diversity and abundance in soils. Proceedings of the National Acad- 
emy of Sciences of the United States of America, 109, 19709-19714. 

Schloss PD, Westcott SL, Ryabin T et al. (2009) Introducing mothur: open-source, 
platform-independent, community-supported software for describing and com- 
paring microbial communities. Applied and Environmental Microbiology, 75, 7537— 
7541. 

Schultz N, JA Silva, Sousa JS et al. (2014) Inoculation of sugarcane with diazotrophic 

bacteria. Revista Brasileira De Ciencia Do Solo, 38, 407-414. 
hapleigh J (2013) Denitrifying Prokaryotes. In: The Prokaryotes (eds Rosenberg E, 


Nn 


Delong E, Lory S, Stackebrandt E, Thompson F), pp. 405-425. Springer, Berlin 
Heidelberg. 
harmin F, Wakelin S, Huygens F, Hargreaves M (2013) Firmicutes dominate the 


N 


bacterial taxa within sugar-cane processing plants. Scientific Reports, 3, 3107. 


e] 


hen T, Stieglmeier M, Dai J, Urich T, Schleper C (2013) Responses of the terrestrial 
ammonia-oxidizing archaeon Ca. Nitrososphaera viennensis and the ammonia- 
oxidizing bacterium Nitrosospira multiformis to nitrification inhibitors. FEMS 
Microbiology Letters, 344, 121-129. 

heneman L, Evans J, Foster JA (2006) Clearcut: a fast implementation of relaxed 
neighbor joining. Bioinformatics, 22, 2823-2824. 

angen JHG, Kerkhoff P (1984) Nitrification inhibitors in agriculture and horticul- 
ture: a literature review. Fertilizer Research, 5, 1-76. 

Soares JR, Cantarella H, Vargas VP, Carmo JB, Martins AA, Sousa RM, Andrade CA 

(2015) Enhanced-efficiency fertilizers in nitrous oxide emissions from urea 


N 


Dn 


applied to sugarcane. Journal of Environmental Quality, 44, 423—430. 

Spang A, Hatzenpichler R, Brochier-Armanet C et al. (2010) Distinct gene set in two 
different lineages of ammonia-oxidizing archaea supports the phylum Thau- 
marchaeota. Trends in Microbiology, 18, 331-340. 

Stieglmeier M, Mooshammer M, Kitzler B, Wanek W, Zechmeister-Boltenstern S, 
Richter A, Schleper C (2014) Aerobic nitrous oxide production through N-nitro- 
sating hybrid formation in ammonia-oxidizing archaea. ISME Journal, 8, 1135- 
1146. 

Stott MB, Dunfield PF, Crowe MA (2009) Class of Chloroflexi-like thermophilic cel- 
lulose degrading bacteria. Google Patents No. CA, 2716728, A1. 

Thauer RK, Jungermann K, Decker K (1977) Energy-conservation in chemotropic 
anaerobic bacteria. Bacteriological Reviews, 41, 100-180. 

Toyoda S, Yano M, Nishimura S-I et al. (2011) Characterization and production 
and consumption processes of N2O emitted from temperate agricultural soils 
determined via isotopomer ratio analysis. Global Biogeochemical Cycles, 25, 
GB2008. 

Urquiaga S, Xavier RP, De Morais RF et al. (2012) Evidence from field nitrogen bal- 
ance and N-15 natural abundance data for the contribution of biological N-2 fixa- 
tion to Brazilian sugarcane varieties. Plant and Soil, 356, 5-21. 

Vandekerckhove TTM, Willems A, Gillis M, Coomans A (2000) Occurrence of novel 
verrucomicrobial species, endosymbiotic and associated with parthenogenesis in 
Xiphinema americanum-group species (Nematoda, Longidoridae). International 
Journal of Systematic and Evolutionary Microbiology, 50, 2197-2205. 


© 2015 The Authors Global Change Biology Bioenergy Published by John Wiley & Sons Ltd., doi: 10.1111/gcbb.12284 


BIOINDICATORS TO IMPROVE C YIELD AND N EFFICIENCY 13 


Varner RK, Keller M, Robertson JR et al. (2003) Experimentally induced root mortal- 
ity increased nitrous oxide emission from tropical forest soils. Geophysical Research 
Letters, 30, 1144. 

Wang Q, Garrity GM, Tiedje JM, Cole JR (2007) Naive Bayesian classifier for rapid 
assignment of rRNA sequences into the new bacterial taxonomy. Applied and 
Environmental Microbiology, 73, 5261-5267. 

Wilson DB (2009) Cellulases and biofuels. Current Opinion in Biotechnology, 20, 295- 
299. 

Zaitsev GM, Tsitko IV, Rainey FA, Trotsenko YA, Uotila JS, Stackebrandt E, Sal- 
kinoja-Salonen MS (1998) New aerobic ammonium-dependent obligately oxalo- 


Supporting Information 


trophic bacteria: description of Ammoniphilus oxalaticus gen. nov., sp. nov. and 
Ammoniphilus oxalivorans gen. nov., sp. nov. International Journal of Systematic 
Bacteriology, 48, 151-163. 

Zhu X, Burger M, Doane TA, Horwath WR (2013) Ammonia oxidation pathways 
and nitrifier denitrification are significant sources of N20 and NO under low oxy- 
gen availability. Proceedings of the National Academy of Sciences of the United States 
of America, 110, 6328-6333. 

Zumft WG (1997) Cell biology and molecular basis of denitrification. Microbiology 
and Molecular Biology Reviews, 61, 533-616. 


Additional Supporting Information may be found in the online version of this article: 


Table S1. Soil physicochemical characteristics of the experimental site. 
Table S2. Chemical composition of the normal vinasse and concentrated vinasse applied as potassium source in soil. 
Table S3. Parameter estimation for fixed effects on global and best models fitted to explain N20 fluxes from soil considering phy- 


sic-chemical parameters and management practices. 


Fig S1 Box-plot representation of ammonium concentration in soil after fertilization. 
Fig S2 Box-plot representation of nitrate concentration in soil after its fertilization. 
Fig S3 Cumulative (%) nitrous oxide emissions from soil under different treatments. 


Fig S4 Relative abundance of Nitrosospira sp. 


Fig S5 Intercept of the negative binomial generalized linear model fitted to verify the effect of treatments on abundance of OTUs. 
Fig S6 Relative abundance of Firmicutes OTUs that significantly explain N20 fluxes. 

Fig S7 Relative abundance of Proteobacteria OTUs that significantly explain N2O fluxes. 

Fig S8 Relative abundance of OTUs belonged to other phyla than Firmicutes and Proteobacteria that significantly explain N,O 


fluxes. 


File S1. Global, nested and best models fitted with microbial taxa to explain N2O fluxes from soil. 


© 2015 The Authors Global Change Biology Bioenergy Published by John Wiley & Sons Ltd., doi: 10.1111/gcbb.12284 


